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Abstract. We introduce the -1 dual Harm polynomials through an appro- 
priate q — y — 1 limit of the dual q-Hahn polynomials. These polynomials are 
orthogonal on a finite set of discrete points on the real axis, but in contrast 
to the classical orthogonal polynomials of the Askey scheme, the -1 dual Hahn 
polynomials do not exhibit the Leonard duality property. Instead, these poly- 
nomials satisfy a 4-th order difference eigenvalue equation and thus possess a 
bispcctrality property. The corresponding generalized Leonard pair consists 
of two matrices A, B each of size N — 1. X :Y f 1. In the eigenbasis where the 
matrix A is diagonal, the matrix B is 3-diagonal; but in the eigenbasis where 
the matrix B is diagonal, the matrix A is 5-diagonal. 



1. Introduction 

Recently new explicit families of " classical" orthogonal polynomials P n (x) were 
introduced [12], [13], [14]. These polynomials satisfy an eigenvalue equation of the 
form 

(1.1) LP n {x) = X n P n (x). 

The operator L is of first order in the derivative operator d x and contains moreover 
the reflection operator R defined by Rf{x) = f(—x); it can be identified as a first 
order operator of Dunkl type written as 

(1.2) L = F(x)(I - R) + G{x)d x R 

with some real rational functions F(x), G{x). The corresponding polynomial eigen- 
solutions P n (x) can be obtained from the big and little q- Jacobi polynomials by an 
appropriate limit q — > — 1. 

In [11] we generalized this approach to the case of Dunkl shift operators. In this 
case the operator L contains the shift operator T + f(x) = f{x + 1) together with 
the reflection operator R: 

(1.3) L = F(x)(I — R) + G{x)(T + R - I) 

(I stands for the identity operator). The rational functions can be recovered from 
the condition that the operator L stabilizes the spaces of polynomials, i.e. it sends 
any polynomial of degree n into a polynomial of the same degree. It can then 
be demonstrated [11] that the polynomial eigensolutions P n {x) of the eigenvalue 
equation (1.1) with the operator (1.3) satisfy the 3-term recurrence relation and 
hence are orthogonal polynomials. In fact, the polynomials P n {x) in this case 
coincide with the Bannai-Ito (BI) polynomials first constructed in [1] (see also 
[101). 
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The BI polynomials thus possess the bispectrality property: they satisfy simul- 
taneously the 3-term recurrence relation (common to all orthogonal polynomials) 

(1.4) Pn+l(x) + b n P n (x) + UnPn-iix) = xP n (x) 

and the eigenvalue equation (1.1). Moreover, in the case when the support of the 
orthogonality measure consists of a finite number of points x s , s = 0, 1, 2, . . . , N, 
the BI polynomials satisfy the Leonard duality property [6], [1], [10]. This means 
that there is a finite difference equation of the form 

(1.5) U s (P n (x s+1 ) - P n (x s )) + V s (P n ( X s — l) Pn{ x s)) — ^nPn{%s) 

with some real coefficients V S ,U S . The difference equation (1.5) is of the second 
order and can be considered as a dual relation with respect to the recurrence rela- 
tion. In fact, the difference equation equation (1.5) is a simple consequence of the 
eigenvalue equation (1.1) with the Dunkl shift operator (1.3) [11]. 

We showed in [11] that the BI polynomials can be obtained by an appropriate 
q — > — 1 limit from the Askey- Wilson polynomials. Correspondingly, the Dunkl 
shift operator L appears in the same limit from the Askey- Wilson difference opera- 
tor [11]. It should be stressed that there are several possibilities in taking the limit 
q — > — 1 of the Askey- Wilson polynomials. Not all of them lead to the BI polynomi- 
als. There is another family of orthogonal polynomials - called the complementary 
BI polynomials in [11] - which can be obtained from the Askey- Wilson polynomi- 
als in the same limit. In contrast to the BI polynomials, the complementary BI 
(CBI) polynomials do not satisfy an eigenvalue equation of the form (1.1). As a 
consequence, the CBI polynomials do not possess the Leonard duality property 
(1.5). 

The main purpose of the present paper is to study the orthogonal polynomials 
which appear in the q — > — 1 limit of the dual q-Hahn polynomials. We call these 
polynomials the dual -1 Hahn polynomials. We derive explicitly basic properties 
and relations for them including a 3-term recurrence relation. The main result is 
the existence of a dual eigenvalue equation (1.1) for the dual -1 Hahn polynomials. 
In contrast to the BI polynomials, the operator L is now 2-nd order with respect to 
the shift operators. As a consequence, the dual -1 Hahn polynomials obey a 5-term 
difference relation on an appropriate grid x s (instead of the 3-term relation (1.5) 
for the BI polynomials). 

2. Dual q-Hahn polynomials 
The dual q-Hahn polynomials [5] 

(2.1) R n {x;a,b,N) = &J q I <?; Q 

depend on 3 parameters a, 6, A, where A is a positive integer. The argument x can 
be parametrized in terms of s as follows 

(2.2) x = q- s +abq s+1 . 

The polynomials R n {x; a, 6, A) satisfy the 3-term recurrence relation 
(2.3) 

A n R n+ i(x; a, 6, N)-(A n +C n -l-abq)R n (x; a, 6, A)+C„i?„_i(x; a, b, N) = xR n (x; a, b, A), 
where 

(2.4) A n = (l-q n - N )(l-aq n+1 ), C n = aq{\ - q n ){b - q^' 1 ). 
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The polynomials R n (x; a, 6, N) are not monic; their monic version P n (x; a, b, N) = 
K n Rn(x;a 7 b, N) — x n + 0(x n ^ 1 ) (with n n an appropriate factor) will obey the 
recurrence relation 

(2.5) P n+1 (x;a,b,N) + b n P n (x; a, b, N) +u n P n -i(x;a,b,N) = xP n (x; a, b, N), 
where 

b n = 1 + abq - A n - C n , u n = A n _iC n . 

The dual q-Hahn polynomials also verify a q-difference equation of second order [5] 
(2.6) 

B(s)R n {x s+1 )+D{s)R n (x s - 1 )-(B{s)+D(s))R n (x s ) = (q- n -l)R n (x s ), s = 0, 1, 2, 
where x s is given by (2.2) and 

(1 - q s - N )(l - aq s+1 )(l - abq s+1 ) 



B(s) = 



(2.7) D(s) = - 



(1 - abq 2 ^ 1 )^ - abq 2s + 2 ) 
aq s ~ N {l - q s )(l - abq s+1+N )(l - bq s ) 



(1 - afeq 2s + 1 )(l - abq 2s ) 

As per [7], the equation (2.6) can be interpreted as a q-difference equation on the 
"q-quadratic grid" x s . 

The dual q-Hahn polynomials satisfy the orthogonality relation 

JV 

(2.8) ^ w s R n {x s )R m {x s ) = K uiux . . . u n S nm , 

s=0 

where the discrete weights are 

(9 q] = (aq,abq,q- N ;q) s 1 - abq 2s+1 JV5 _ s ( s _i)/ 2 

{ ' S (q,abq»+ 2 ,bq;q) s (l-abq)(-aqy q 

and the normalization constant is 

(abq 2 ;q) N N 

(bq;q)N 

We used the standard notation for the q-shifted factorials [5] . 

When q — > 1 and a = q a , b = q^ the dual q-Hahn polynomials become the 
ordinary dual Hahn polynomials [5] 

/— n — s s + 1 + a + 
(2.10) W n {x s -a,b 1 N) = z F 2 [ ' T tr I 1 



a + l,-N 
where 

x s = s(s + a + + 1). 
These polynomials satisfy the three-term recurrence relation 

(2.11) A n W n+1 (x s ) - (A n + C n )W n (x s ) + CnWn-iixs) = x s W n (x s ), 
where 

A n = (n - N)(n + a + 1), C n = n(n - - N - 1). 
The corresponding monic dual Hahn polynomials W n (x) obey 

(2.12) W n+ i(x s ) + b n W n (x s ) + u n W n -i(x s ) = x s W n (x s ), 
where 

(2.13) u n = An-iCn = n{n-0-N- l)(n — N — l)(n + a), b n = -A n - C n . 



4 



SATOSHI TSUJIMOTO, LUC VINET, AND ALEXEI ZHEDANOV 



3. A LIMIT q^f — 1 OF THE RECURRENCE RELATIONS AND THE ORTHOGONALITY 

PROPERTY 

We wish to consider a limit q — > — 1. We will assume that a — > ±1 and b — > 
±1 when q — V — 1. We want to obtain a nondegenerate limit of the coefficients 
A n {l + q)^ 1 , C„(l + q)^ 1 for <j — > —1. This means that both these limit coefficients 
should exist and be nonzero for all admissible values n = 1, 2, . . . , N. It is hence 
easily seen that necessarily, we must have ab — > I. Two situations have to be 
considered separately: 

(i) when N — 2,4,6,... is even, the nontrivial q — > — I limit then exists iff 
a — >• 1, & — >• 1. It is then natural to take the parametrization 

(3.1) q = -e e , a = e- aE , b = e- fjE , e^O 

with real parameters a,/3. Dividing the recurrence relation (2.3) by q+1 and taking 
the limit £->0we obtain the recurrence relation 

(3.2) (v.; «. A + ^"'^i-V (y,; «. A 

(4" 1 ) + C^- 1 ))^- 1 )^; a, /3, JV)) = VsR^ {y s ;a, /3, TV), 
where the grid y s has the following expression 
, x J —a — (3 + 2s + I if s even, 

{3 - 6) Vs ~\ a + /3-2s-l if « odd 

or, equivalently, 

(3.4) y s = (-l) s {I- a- /3 + 2s), s = 0, 1, . . . , N. 

The recurrence coefficients are 
(3.5) 

^(-i) f 2(n — N) if n even ^(-i) f ^ n even 



2(ra+l-a) if n odd' " [2(N+l-/3-n) if n odd 
The corresponding monic —1 dual Hahn polynomials satisfy relation (2.5), where 

(3 6) u(-V - A^tf-V ~ ! 4n(«-«) if « even 



4(A r -n+l)(n + /3-N-l) if n odd 
and 

(37) b (-V-i- a -B-A(-V-C(-V-l 2N + 1 ~ a ~ ? if n cvcn 
K - 1 a P An °" "\-27V-3 + a + /3 if n odd 

It is convenient to introduce the "^-numbers" 

(3.8) [n]^ = n + M (l-(-l) n ), 

which appear naturally in problems connected with the Dunkl operators [8]. One 
can then present the recurrence coefficients in the compact form 

(3.9) u^-V =4[n]t[N-n + l] v , = 2([n] ( + [N - n] v ) + I - a - /3, 
where; 

, 0-N-1 a-N -1 

?= 2 = 2 ' 

It is seen that u = u^+i = as required for finite orthogonal polynomials. The 
positivity condition u n > 0, n = 1, 2, . . . , N is equivalent to the conditions 

(3.10) a > N, > N. 
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These polynomials arc orthogonal on the finite set of points y s 

N 

(3.ii) Y, WsR n 1 Hys)R ( m 1 Hys)^n u[- 1) ui- 1) ...ui- 1 H nm , 

s=0 

where the discrete weights are defined as 
(3.12) 

HV/2) fl (l-a/2) 8 (l-a/2-/3/2) fl N 
W2s ~ { l) si (l-/3/2) s (7V/2 + l-a/2-/3/2) s ' s - u >^--- 2 

and 
(3.13) 

_ (-7V/2) S+1 (l- a /2) a (l- a /2-/3/2), iV 

25+1 " ( j s! (1-/3/2)^/2 + 1- a/2 -/?/2) s+1 ' s - U ' i '--- 2 ^ 
The normalization coefficient is 

d-S±S) 

1 ) N/2 



(3.14) k 



'-«). 



(i 

V "/N/2 

Assume that a = TV + ei, /3 = N + e 2l where e lj2 are arbitrary positive pa- 
rameters. This parametrization corresponds to the positive condition for the dual 
-1 Hahn polynomials. Then it is easily verified that all the weights are positive 
w s > 0, s = 0,1,..., N. 

Moreover, the spectral points y s are divided into two non-overlapping discrete 
sets of the real line: 

{1 - 5, -3 - 5, -7 - 5, ... , -2N +1-5} 

and 

{1 + 5,5 + 5,9 + 5,..., 2N -3 + 5}, 
where 5 = e\ + e 2 > 0. The first set corresponds to y s with even s and contains 
1 + N/2 points; the second set corresponds to y s with odd s and contains N/2 
points. 

(ii) when N = 1, 3, 5, . . . is odd, a nontrivial q -+ — 1 limit also exists iff a -+ 
— 1, b — > — 1. We take the parametrization 

(3.15) q = -e e , a = -e a£ , b = -e' is , e^O 

with real parameters a, /3. Dividing again the recurrence relation (2.3) by q+1 and 
taking the limit e — > 0, we obtain the recurrence relation (3.2) where the grid y s is 
defined as 

. „ ,-]3 + 2s+l if s even 

and the recurrence coefficients given by 
(3.17) 

— 2n if n even 
2(/3 + iV-n+l) if n odd ' 

The corresponding monic —1 dual Hahn polynomials satisfy the standard relation 
(2.5), where 

, (-1)^-1) _ / 4n(N+l-n) if n even 

■ + n){p + N + l-ri) if n odd 



( a + 
\-a 
cients giv 

(2(a + n + l) if n even f 
\ 2(n-JV) if n odd ' " \: 

The corresponding monic —1 dual Hahn 
(2.5), where 

(3-18) tt ( r i) = ^-i)q(-i) = | 4(a _ 
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and 

- 1 — a + f3 if n even 
— 1 + a — ft if n odd 



(3.19) b(-V = l + a + f} - - = | 



Again, as in the case of even N, it is possible to present the recurrence coefficients 
in the compact form 

(3.20) ui-V =4[nk[JV-n+l]„ = 2([n} 6 + [N - n} n ) - 2N - 1 - a - f3, 

with £ = a/2, r\ = j3/2. It is seen that in both cases: N even and N odd, the 
recurrence coefficients of the dual -1 Hahn polynomials are presented in the unified 
form (3.9) or (3.20) with the difference only residing with the parameters £,77. 

It is seen that uq = u^+i = as required for finite orthogonal polynomials. The 
positivity condition u n > 0, n = 1, 2, . . . , N is equivalent either to condition 

(3.21) a > -1, P > -1 

or to condition a < —N, (3 < —N. In what follows we shall use only condition 

(3.21) . 

The polynomials i?l ^(2;) are orthogonal on the finite set of points y s 

N 

(3.22) J2 w - R n 1 Hvs)R^ 1 Hys) = ^u { 1 - 1) ui- 1) ...ui- 1 U nm , 

s=0 

where the discrete weights are defined as 
(3.23) 

' _ ( , )S (-(N-1)/2) S (l/2 + a/2) s (l + a/2 + p/2) s , _ ni9 * 

2s [ ' s\ (l/2 + /V2) s (7V/2 + 3/2 + a /2 + /V2y U ' i ' / '"'-' 

and 
(3.24) 

_ (-(N-l)/2) s (l/2 + a/2) s+1 (l + a/2 + p/2) s s _ nl9 

2S+1 [ ' s\ (l/2 + /3/2) s+1 (7V/2 + 3/2 + a /2 + /3/2) s ' U ' i,A -' 

The normalization coefficient is 
(3.25) « = V ){N+1) ' 2 



(JV+l)/2 



Assume that a = —1 + ei, /3 = —1 + £2, where ei,2 are arbitrary positive 
parameters. This parametrization corresponds to the positive condition for the 
dual -1 Hahn polynomials for TV odd. Then it is easily verified that the weights are 
positive w s > 0, s = 0, 1, . . . , N. 

Moreover, the spectral points y s are divided into two non-overlapped discrete 
sets of the real line: 

{-1 - (5, -5 - (5, -9 - (5, ... , -2N + 1-5} 

and 

{-1 + 5,3 + 5,7 + 5,..., 2N -3 + 5}, 
where 5 = t\ + €2 > 0. Both sets contain (N — l)/2 points. 
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4. Explicit expression in terms of the ordinary dual Hahn 

polynomials 

In this section, we derive an explicit expression for the -1 dual Hahn polynomials 
in terms of the ordinary dual Hahn polynomials. 

Consider first the case of even N — 2,4,6, Introduce the "shifted" monic 

-1 dual Hahn polynomials R n {x) = Rn X \x — 1) From formulas (3.6), (3.7), we 
conclude that these polynomials satisfy the recurrence relation 

(4.1) Rn+lix) + {-l) n TR n (x) + U n Rn-l(x) = xR n (x), 

where 

T = 2N + 2-a~(3 

and u n are given by (3.6). 

A recurrence relation of the type (4.1) leads to orthogonal polynomials R n (x) 
which are very close to symmetric orthogonal polynomials. Using methods devel- 
oped in [2] and [13], we can introduce a pair of monic orthogonal polynomials P n (x) 
and Q n (x) by the formulas : 

(4.2) R 2n {x) = P n (x 2 ), R 2n+1 (x) = {x- T)Q n (x 2 ). 

It can easily be shown that the polynomials P n (x) and Q n (x) satisfy the following 
recurrence relations (it is assumed that uq — 0) 

(4.3) P n +\{x) + (u 2n + U2n+i + r 2 )P n (x) + u 2n u 2n -iPn-i{x) = xP n (x) 
and 

(4.4) Q n +i(x) + (u 2n+2 + u 2n+1 + T 2 )Q n (x) + M 2 r l U2n+iQr l -i(a;) = xQ n (x), 
and moreover that the polynomials are connected by the Christoffcl transform 

(4.5) Q n (x) = P " +1(X)+M2 " +1P " (X) . 

X — T A 

It is also easily seen that both P n (x) and Q n {x) are ordinary dual Hahn polynomials. 
We hence have the following explicit expression: 



(4.6) R 2 l 1 \x-l)=^hF 2 [^ n 'V, f,,7 a 1 |l), n = 0,l,2,... 

and 



_N_ i _ a 
' 2 ' 1 2 



(4.7) R^%(x - 1) = 7 «0r - r) 3 F 2 ^ „ 1 I lj , « = 0,1,2,..., 

where 77 = 1/2 — (a + /3)/4 and the normalization coefficients are 

7 i 0) - 16"(-iV/2)„(l - a/2)„, ^ - 16"(1 - N/2) n {l - a/2) n . 
Quite similarly, for the odd N = 1, 3, 5, . . . we find 

(4.8) R 2 - 1 \x-l)=^hF 2 (^ n ^^ 1 ^ \lj, n = 0,l,2,... 
and 

(4.9) Rt 1 \x-l)=^\x + a-P),F 2 (^\ + J^~ i \l\, n = 0,l,2,.... 
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where 

, . «±£±», 7 » . ur (i^)^ (2±i)_, - 16- (2±2 

Some of these polynomials have appeared in [9] , [4] in the context of quantum spin 
chains. 

5. Difference equation 

Consider the following operator L defined on the space of functions f(s) that 
depend on a discrete variable s: 

(5.1) Lf(s) = B(s)f(s + 1) + D(s)f(s - 1) - (B(s) + D(s))f(s), s = 0, 1, 2, . . . , 

where B(s),D(s) are given in (2.7). Manifestly, the difference equation (2.6) means 
that the dual g-Hahn polynomials are eigenfunctions of the operator L 

(5.2) LR n (x(s)) = \ n R n (x( s )) 
with the eigenvalues 

(5.3) A„ - q~ n - 1. 

When q — > 1 we obtain the difference eigenvalue equation for the ordinary dual 
Hahn polynomials 

(5.4) LiWnixis)) = -nW n (x(s)), 

where the operator L\ can be obtained from L as L\ = lim 9 _>i L(q— Explicitly 
[5] 

(5.5) U = B x (s)f{s + 1) + Di (*)/(« - 1) - + 
where 

n = (s + a + /3 + l)(s + a + l)(jV- S ) = s(s + ^)(a + Q + /3 + iV+l) 

l(S) ~ {2s + a + p + l)(2s + a + /3 + 2) ' l(Sj ~ (2s + a + + l)(2s + a + 0) ' 

When we try to perform a similar procedure for the limit q — > — 1 we encounter a 
problem. Indeed, it is easily seen that L(l + ?) _1 does not have a nondegenerate 
limit as q — > — 1 . We thus cannot obtain an eigenvalue equation in 3-diagonal form 
like (5.4) for the dual -1 Hahn polynomials. 

Nevertheless we observe that the operator (L 2 + 2L)(1 + q)^ 1 does survive in 
the limit q — > — 1. We hence have the following eigenvalue equation for the dual -1 
Hahn polynomials 

(5.6) HR(- 1 \y s ) = 2nKi- 1 \y s ) 1 
where the grid y s is defined by (3.3) or (3.16) and 

H = lim (L 2 + 2L)(l + q)- 1 . 

3->-l 

The operator H obviously is 5-diagonal, i.e. 
(5.7) 

^/( s ) = C/ 2 ( s )(/( s +2)-/( s ))+C/ 1 ( s )(/( S +l)-/( S ))+y 2 (s)(/( S -2)-/( S ))+y 1 ( S )(/(s-l)-/( S )) 
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The explicit expressions for the coefficients Ui(s), Vi(s) depend on the parity of N. 
For even N = 2, 4, 6, . . . they are: 



(5.8) U 2 (s) 



(5.9) Ui(a) 



(5.10) V 2 (s) 



o (a-s-2)(/3+a-s-2)(N-s) -r 

Z (q+^-2s-2)(q+,3-2s-4) 11 6 eVt3Ii 

O (g- s -l)(ff + a - a -l)(JV- s -l) - ,, 

Z (a+/9-2s-2)(a+,3-2s-4) 11 6 OUU 



1 ( Q+/ 3_2 s )( Q+( 3-2s-2)(a+/3-2 S -4) 11 S eVen 
, ( a -s-l)03+ a -s-l)(2iV+2-a-,8) . f 
* (a+0-2s)(a+/3-2s-2)(a+0-2s-4) 11 * UUU 



^ (ct+/ 3_ 2s +2)(Q+/3-2 ;S ) 11 S eVen 
2 ( Q +/3-2 S +2)(a+0-2 S ) 11 S 0dCl 



-s)(2 N+2-a-£) 



(5.11) Vi(s) = 



(a+l3-2s)(a+P-2s-2)(a+l3-2s+2) ^ S eVen 
(j3+a)(a-/3)(-a-i3+iV+s+l) j f g ^ 



(a+^-2s)(a+,3-2s-2)(a+ / 3-2s+2) 

and for odd N = 1, 3, 5, . . . 



(5.12) t/ 2 ( S ) = 



(5.13) E/i(a) = 



(5.14) y 2 ( a ) 



(5.15) 14 (s) 



O ( a +/3+ S +2)(a+s+l)(iV-s-l) jf 

_Z (a+/3+2s+2)(a+,8+2s+4) 11 S eVen 

r, (a+s+2)(a+/3+ S +l)(jV- S ) . f 

~ Z { a +fi+2 s+2)( a +f3+2s+4) 11 S 0dd 



o (a+/3)( a +s+l)( a +/3+2iV+2) . f 

Z («+^+2s)( Ct +/3+2 S +2)(a+ / 3+2 ;S +4) 11 6 eVW1 



(a-P)(N-s)(a+P+s+l) 
(a+P+2 s) (a+13+2 s+2)(a+/3+2 s+4) 



if s odd 



2 s(P+s-l)( a +P+N+s+l) . f 

Z (a+/3+2s-2)(Q+ff+2s) 11 6 eVe11 

Z (a+£+2s-2)(a+/3+2s) 11 6 OUU 



* (a+/3+2s)(a+/3+2 S +2)(a+/3+2 S -2) 11 6 eVt311 

(/3+ S )(g+ff)(g+ff+2jV+2) . f 
Z (a+/^2 S )(a+/3+2s+2)(a+/3+2 S -2) 11 S 



We thus have a difference equation for the dual -1 Hahn polynomials in the form 

u 2 (s) (i^W) - Ri-'Hvs)) + Ui(s) (Ri-'Hvs+i) Rk-'Kvs)) + 

V 2 (s) (4- 1) (y,- 2 ) - R { - 1} (ys)) + Vi(s) (fl^Gfe-i) - ^ InR^iys), 

where the coefficients Ui(s), Vi(s) are provided in the above formulas. 

6. Another form of the difference equation 

The difference equation for the dual -1 Hahn polynomials can be presented in a 
more compact form if one notices that the grid y s satisfies the relations 

, C1 , / ~Vs T 2 if s even 

y —y s ± 2 if s odd 
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This property implies that the difference equation can be written as 

E 1 (x) (r^Hx + 4) - R^Hx)) + E 2 (x) (R^Hx 4) - R^Hx)) + 

GiW (R^K-x 2) - Rt'Kx)) + G 2 {x) (Rir^-x + 2) - R^Hx)) = 2nR^\x) 
or, in operator form as 

(6.2) HR[- 1 Xx) = 2nR< n - 1 \x). 

The operator H in (6.2) reads 
(6.3) 

H = E 1 (x)T i +E 2 (a;)T- 4 +Gi (x)T 2 R+G 2 (x)T- 2 R- {Ei (x)+E 2 (x)+G 1 (x)+G 2 {x))I, 

where the operators T and T~ x are the standard shift operators: T j f(x) = f(x + 
i)> 3 = 0, ±1, ±2, . . . , and R is the reflection operator Rf(x) = f(—x). I denotes 
the identity operator. 

The functions Ei(x), Gi(x), i = 1, 2 are simple rational functions in x. 

For even N — 2, 4, . . . we have 

(x + 3- a + /3) (x + 3 - a- (i) {x -1 -2N + a + /3) 



£ 2 (x) = - 



4(x + l)(3 + ar) 
(a-l + ^ + a;)(a-l-/3 + a;)(ar-l + 2iV-a- / 8) 



Gi(x) 
G 2 (x) = 



4 (a; - 1) (a; - 3) 
_ (a 2 - i3 2 ){x + a + P-2N -1) 
~ (x 2 - 1) (x + 3) ' 

(2JV + 2 - a - p) ({x + a - l) 2 - f3 2 



(x 2 -l)(x- 3) 
For odd N = 1,3,5,... 

{x + a + f3 + , i)(x + a-f3 + l){x-a-l3-2N + l) 



Ei(x) = 



E 2 (x) = - 



4(x + l)(x + 3) 
(x-a-(3-l)(x + P-a-?,)(x + a + P + 2N+l) 



G 2 (x) = 



4 (x - 1) (a; - 3) 
_ (g + g) (g + g + 2 + 2 N) (x + 1 + a - p) 

(l-a; 2 )(x + 3) 
(a - P) {x - a - 13 - 1) (x + a + /3 + 2 N + 1) 



(l-x 2 )(x-'S) 

The form (6.2) of the difference equation is preferable because we here have the 
operator H acting directly on the argument of the polynomials. The operator H 
belongs to the class of Dunkl shift operators: it contains both simple shifts and 
the reflection operator R. Moreover, the operator H preserves the space of poly- 
nomials: it transforms any polynomials of degree n into a polynomial of the same 
degree n. Operators of this kind were considered in the theory of Bannai-Ito poly- 
nomials [11]. However, in contrast to the Bannai-Ito situation, the operator H in 
the present case is of second order, while the Bannai-Ito polynomials are eigenfunc- 
tions of a Dunkl-shift operator of the first order [11]. Equivalcntly, this means that 
for generic values of the parameters a,f3, the dual -1 Hahn polynomials are eigen- 
vectors of a 5-diagonal matrix (while the Bannai-Ito polynomials are eigenvectors 
of a 3-diagonal matrix). 



DUAL -1 HAHN POLYNOMIALS 
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